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Abstract 



Finite temperature properties of a non-Fermi liquid system is one of the 
most challenging probelms in current understanding of strongly correlated 
electron systems. The paradigmatic arena for studying non-Fermi liquids is 
in one dimension, where the concept of a Luttinger liquid has arisen. The 
existence of a critical point at zero temperature in one dimensional systems, 
and the fact that experiments are all undertaken at finite temperature, implies 
a need for these one dimensional systems to be examined at finite tempera- 
ture. Accordingly, we extended the well-known bosonization method of one 
dimensional electron systems to finite temperatures. We have used this new 
bosonization method to calculate finite temperature asymptotic correlation 
functions for linear fermions, the Tomonaga-Luttinger model, and the Hub- 
bard model. 
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I. INTRODUCTION 



The many-body problem in condensed matter physics generally relies on the machinery 
of perturbative quantum field theory to obtain solutions to model systems. Solutions to 
the field theories are obtained in terms of Green's functions, or equivalently, correlation 
functions. The Kubo formulae relate the retarded Green's functions of a model to the linear 
response of the system due to external fields. The experimental responses of a real system 
may then be compared to the solutions for the model. 

Experiments on many-body electron systems normally measure properties at energy 
scales small compared to the Fermi energy. This means that only a few degrees of free- 
dom in the system are excited, and only the low energy sector of a model need be compared 
against experiment. 

The low energy, long distance physics also determines any long range order and coop- 
erative phenomena present in a system. Divergences in certain correlation functions are 
indicative of the presence of ordering phase transitions in the system, such as ferromag- 
netic transitions or the Cooper pairing superconductivity phase transition. The field theory 
methodology has been extremely successful, with the properties of a number of simple met- 
als accurately modeled by using Landau's Fermi liquid theory. However, many materials 
have emerged that have non-Fermi liquid properties. 

The paradigmatic arena for the study of non- Fermi liquid behaviour is in one dimensional 
many body systems. Simple dimensional arguments show that whilst finite momentum 
transfer interactions become irrelevant at low energies in higher dimensions leading to Fermi 
liquid behaviour, in one dimension the finite momentum transfer interactions are marginally 
relevant at all energies. 

Normal perturbative methods of solution can no longer be justified in one dimension. 
This is because the non-interacting one dimensional fermion gas is unstable against the 
switching on of interactions. This is a form of the orthogonality catastrophe, where the 
interacting ground state is orthogonal to the non-interacting ground state. For this reason, 
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non-perturbative methods of solution are now being examined. 

Ground state properties of one- dimensional Fermi liquids are very special in that they 
retain a Fermi surface (if defined as the set of points where the momentum distribution or its 
derivatives have singularities) enclosing the same /c-space volume as that of free fermions, 
in agreement with Luttinger's theorem. However, there are no fermionic quasi-particles, 
and their elementary excitations are instead bosonic collective charge and spin fluctuations 
dispersing with different velocities. An incoming electron decays into such charge and spin 
excitations which then spatially separate with time (charge-spin separation). 

Spin and charge separation is seen in all one-dimensional electron systems and it seems 
that it is not related to other characteristics of the one-dimensional phenomena. The reason 
for this is simple: spin and charge separation is due to scattering at the same Fermi point. 
Accordingly, it will not give any contributions to effects which are due to scattering be- 
tween two different Fermi points, e.g., the renormalized correlation function exponents. The 
correlations between these excitations are anomalous and show up as interaction-dependent 
nonuniversal power-laws in many physical quantities where those of ordinary metals are 
characterized by universal (interaction independent) powers. 

These properties are generic for one-dimensional Fermi liquids but particularly prominent 
in a one dimensional model of interacting fermions proposed by Tomonaga (Tomonaga 1950), 
Luttinger (Luttinger 1963) and solved exactly by Mattis and Lieb (Mattis and Lieb 1965). 
All correlation functions of the Luttinger model can be computed exactly, so that one has 
direct access to all physical properties of interest. The notion of a Luttinger liquid was given 
by Haldane (Haldane 1981) to describe the universal low-energy properties of gapless one 
dimensional quantum systems, emphasizing that an asymptotic {uj — > 0, q — > 0) description 
can be based on the Luttinger model in much the same way as the Fermi liquid theory 
in three dimension is based on the free Fermi gas. The basic ideas and procedures had 
been discussed earlier by Efetov and Larkin (Efetov and Larkin 1975) but passed largely 
unnoticed. 
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A. Thermodynamic Bethe Ansatz 



The first approach to extend Bethe Ansazt to finite temperatures was developed by 
Yang and Yang (Yang and Yang 1969; Yang 1970) for the one dimensional boson gas with 
5-function interaction. The method is based on the analysis of the spin wave spectrum of the 
Bethe Ansatz equations. The first analysis of this spectrum took into account only real roots 
of the Bethe Ansatz equations (Cloiseaux and Pearson 1962; Griffiths 1964). But in general 
the roots are complex, as first observed by Katsura (Katsura 1965; Ovchinnikov 1969), 
which give bound states and may be grouped in various "strings". This "string" method 
was successfully implemented in all cases where the ground state contains bound pairs, such 
as: i) the |A| = 1 Heisenberg spin-chain (Takahashi 1971a). it) the |A| > 1 Heisenberg-Ising 
(XXZ) spin-chain (Gaudin 1971; Takahashi and Suzuki 1972) (extended also for general XYZ 
model (Takahashi and Suzuki 1972)). Hi) the ID attractive 5-function interacting electron 
systems (Takahashi 1971b; Schlottmann 1993). Hence also for the U < Hubbard model 
and half-filled U > Hubbard model (Takahashi 1972); and iv) models with effective one 
dimensional attractive ^-function interactions such as the Kondo model (Filyov, Tsvelick and 
Wiegmann 1981; Andrei and Lowenstein 1982; Tsvelick and Wiegmann 1983a), the Coqblin- 
Schrieffer model (Tsvelick and Wiegmann 1982a; Rasul 1982; Tsvelick and Wiegmann 1983a; 
Andrei, Furuya and Lowenstein 1983) and the nondegenerate Anderson model with infinite 
U (Filyov, Tsvelick and Wiegmann 1982; Tsvelick and Wiegmann 1982b; Kawakami and 
Okiji 1982; Tsvelick and Wiegmann 1983a; 1983b, 1983c; Okiji and N. Kawakami 1983; 
Schlottmann 1984). 

The string (or Thirring) hypothesis is valid for models exhibiting Luttinger liquid ground 
state properties, repulsive 5-function interacting models (Takahashi 1971b, 1972; Lai, 1971, 
1973), and the partially filled U > Hubbard model (Takahashi 1972). As such, the same 
properties apply. The Bethe Ansatz equations of the Hubbard model have three types of 
roots: i) real k, ii) k and A string solutions and in) pure A strings. The existance of the real 
k roots shows that even at finite temperature pure holon type excitations, and hence spin 
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- charge separation, exist. For more details, see Korepin, Bogoliubov and Izergin (1993). 
However, even in cases where the Bethe Ansatz method is applicable, correlations functions 
cannot be calculated and other methods, i.e. conformal field theory or bosonization has to 
be implemented. The later is our main focus. 

B. Finite Temperature Formalism 

The static properties of an equilibrium system at given temperature can be ascertained 
by deriving the partition function of the system, giving standard thermodynamic properties, 
(Garrod 1995). To calculate the dynamical properties of systems at finite temperature we 
must adopt a formalism that includes the statistical mechanics of a system for a given 
temperature. 

The finite temperature expectation values of a field theory are no longer ground state 
expectations, as they are for theories at T = 0, but expectations taken over an ensemble of 
states. The operator corresponding to the observable must be averaged over an ensemble of 
basis states, with the weighting of each state determined by its energy. Mathematically this 
is expressed by, 

(A) p = ^Tr e^A , (1) 

where Tr is the trace over the basis states, and Z is the partition function Z = Tr e~^ H . 

A number of finite temperature formalisms have been developed, mainly relying on a 
mathematical correspondence between quantum field theories in D dimensions, and the 
statistical mechanics of systems in D+l dimensions. 

The simplest finite temperature formalism is the imaginary time, or Matsubara formal- 
ism. In this formalism the equilibrium statistical mechanics of the system is equivalent to an 
imaginary time evolution. The solutions can then be determined by standard methods using 
an "imaginary time" variable r = it, where r is an inverse temperature. Unfortunately, this 
leads to the loss of any information about the dynamics of the system, as the time variable 
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has been replaced by the temperature. Therefore, only static equilibrium properties of the 
system can be analysed. Dynamic properties, such as spin-charge separation, are no longer 
observable in this formalism. 

Dynamic finite temperature correlation functions may be calculated from the results 
of the Matsubara method, however, by using the spectral density function (Inkson 1984). 
The "real time" correlation functions are determined by analytically continuing the spectral 
density function. From the Matsubara method the spectral density function is known at 
infinitely many points on the imaginary line, the Matsubara frequencies iu n . If the complete 
spectral density function can be obtained, then the real time correlation functions can be 
given for arbitrary time and temperature. 

Using the bosonization technique, the trace over a basis may be done using a basis of 
boson number states. A direct evaluation of the expectation values may thus be obtained 
easily, without the use of Matsubara frequencies and analytic continuation. Properties such 
as spin-charge separation may then easily be studied at finite temperature. 

Hereafter, we review the bosonization of certain one dimensional models, and calculate 
the correlation functions, and some Green's functions, for these one dimensional models at 
finite temperature. 

II. BOSONIZATION 

Bosonization of a quantum field theory describes a method by which fermionic operators 
in the theory, obeying anti-commutation relations, are replaced by bosonic operators obeying 
commutation relations. 

Replacement of one field theory with another would appear to merely replace the problem 
of solving the original fermionic field theory with the problem of solving a bosonic field 
theory. The usefulness of this technique was realised when it was discovered that certain one 
dimensional interacting fermionic field theories, were equivalent to non-interacting bosonic 
field theories. While the solution to interacting field theories is often difficult, and generally 
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only perturbative methods are available to obtain a solution, the solution to a non-interacting 
field theory is well known and exact (Abrikosov, Gorkov and Dzyaloshinski 1963; Inkson 
1984). 

Given the exact solution to the bosonic field theory, the properties of the fermionic theory 
can be calculated by use of a "bosonization dictionary". For example, correlation functions 
in terms of fermionic operators can be re-written as expectations for bosonic operators, and 
solved. 

Bosonization may be derived by defining a set of boson fields (f>(x), and their conjugate 
momenta fields H(x), with given commutation relations, and then determining the com- 
mutation relations and Green's functions of the exponentials of these fields. The Green's 
functions and commutation relations are fermionic in nature, and an identity between the 
boson fields, <f>(x) and IL(x), and fermion fields ip(x) is made (Shankar 1995; Gulacsi 1997). 

Field theoretical bosonization evolved out of solutions to models in the area of high 
energy physics. Indeed, many of the models used in the particle physics arena have a direct 
correspondence to models used in condensed matter physics. Some relevant examples of this 
correspondence shall be pointed out later on. 

A. Bosonization Via Fock Space Identification 

Constructive bosonization, in terms of a Fock space identification of operators, gives a 
simpler, more transparent derivation of the boson representation of fermion operators. An 
explicit picture of the physics corresponding to the boson representation may be obtained 
directly from the formalism. We adopt this procedure hereafter, not just because of its 
simplicity, but mainly because it allows a straightforward extension to finite temperature. 

In this section we derive the Fock space identity between fermions and boson coherent 
states, closely following the pedagogical article of Schonhammer and Meden (1996); Delft 
and Scholler (1998). From this derivation we obtain the bosonization dictionary we require 
for the solution of the models included in later chapters. 
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1. Requirements for Bosonization 

The fermionic system must be expressible by a countable set of operators, with un- 
bounded spectrum, that obey fermionic anti-commutation relations. That is, operators Cfc jCr 
with k G [—00, 00], a a species identifying label, and, 

{<w,4<t} = <w<W' • (2) 

The species identifying label can be just the electron spin a = {f, !}, or in the Luttinger 
model the spin combined with the left moving L and right moving R species labels. 

The treatment of fermions on different sections of the Fermi surface as different species 
can be used to generalize this method to higher dimensions. The success of this approach 
in one dimension hinges on the assumption that there is no mixing of fermion species, due 
to the large 2kp jump in momentum between the points of the Fermi surface. 

At this point there are no requirements on the type of Hamiltonian of the system, as the 
operator identity is a property of the Fock space. 

2. The Vacuum State and Normal Ordering 

In order to regularize all the expressions of operators within the theory, a zero particle 
vacuum state |0) for the fermions is introduced. This vacuum state is then used to define 
normal orderings of all the operators. The vacuum state is defined by, 

Cfc<r|0) = ; for k > , 
4|0) = ; for k < . 

Hence, for any combination of the fermion operators we may define the normal ordering 
of such operators by, : ABC : = ABC - (0\ABC\0). 
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3. Defining the Boson Operators 



The boson operators b qa and b qa are defined on the Fock space as linear combinations of 
particle-hole excitations, 

/2~ 00 



j ~K+q,a~K<j 1 ~qa W r C k-q,cr ( 

" k= — oo V y fc= — 00 



fet _ ./f^. V r f r ft 

<?°" \/ T^Y ^ ^k+q,a^ka 5 °V 



fc<T ) 



(3) 



where the momentum g = ^ for n a positive integer. The boson operators are propor- 
tional to the q 7^ Fourier components of the fermion density , ipj T (x)'^ a (x), with the q = 
component related to the fermion number operator N. 
The operators obey the bosonic commutation relations, 



(4) 



where the last equality requires normal ordering and the existence of the negative k states, 



^k+q-q',cr^k,<j ^k+q,a^k+q' ,a 



= f q Saa'S qq/ Ek [(0\cl a C ka \0) - (0\c{ +qta C k+q JO) 



4- Partitioning the Fock Space 

The fermionic Fock space can be interpreted as a direct sum of subspaces, each of defi- 
nite particle number N. Excitations within these subspaces are then bosonic particle-hole 
excitations, with the b g and b q operators acting as annihilation and creation operators, re- 
spectively. As first noted by Overhauser (1965), the bosonic particle-hole excitations form 
a complete basis of states for each N particle subspace. 

The operators that connect the subspaces together are the unitary Klein operators F t 
and F, which increase and decrease the total particle number by one, respectively. In this 
way, the Klein operators act as raising and lowering operators with respect to the fermion 
number operator. They obey the consequent commutation relations, 

[N a , Fl] = 5 aa ,Fl, , [N 9 , F a ,} = S^F, . (5) 
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The operators were first constructed in terms of the bare fermions by Haldane (1979), 
and an overview of their properties given by Haldane (1981); Delft and Scholler (1998). For 
the purposes of this paper we simply describe remaining commutation relations obeyed by 
the Klein operators, which are, 

{Ft, F a ,} = 2<W , {H, fy) = it, K>} = , 

[St,, F%\ = [b qt7 , F%\ = p qa , F a ,\ = [b qa , F a ,\ = . (6) 
We also note that because the Klein operators are unitary, we have, F^F a = F a F^ = 1. 

5. Fermion Operators as Coherent States 

The fermion field operators ip a (x) and ipl(x) are defined in terms of the Fourier compo- 
nents of the c ka and c\ a operators, respectively, 

/ C\ OO IQtt 00 

' fc=— oo ' k=— oo 

The fundamental relation that makes bosonization work is that the state produced by 
the fermion operator ip a (x) on an N particle ground state | AT ) is an eigenstate of the boson 
operator b q a . This may be derived from the commutation relations, 

[b qa ,^)\ = -^>\j Tq e iqx ^Ax) , K^A*)\ = -^^—e-^A*) , (8) 

and the fact that any N particle ground state |JV ) vanishes when operated on by 6 ?(7 . This 
allows the representation of the fermion operator ip a (x) acting on any N particle ground 
state as a coherent state representation of boson operators. Thus, 

2tF .., 



j a (x)\N ) = F a X a (x) exp ( - J2 \J^% j \N ) , (9) 



q>0 



where the Klein operator F a is required to reduce the particle number, and X a (x) is a phase 
operator, 

Xa ( x ) = ^Ee-iW-fr)* . (10 ) 
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The boson field operators ip a (x) and (fi^ix), are introduced to simplify the notation, and 
are simply, 

/ O-jr- . IOtt 

M*) = - £ \ J-/ qXh ^ > ¥>£(*) = " E \ -Tn e ~ i9Xb l ■ (") 
4>0 V ^ g>0 V ^9 

and allow us to write (a?) |iV ) = F CT A< T (:r)e ¥V ( a: )|iVo)- 

6. The Boson Representation of ip a (x) 

Bosonization is only useful because properties of the fermionic system can be described 
in terms of certain boson fields, whose properties are simpler to calculate than those of 
the original fermion fields. The correlation functions, and hence Green's functions, of the 
fermionic system contain combinations of single particle fermion operators, hence we need 
to express these operators in terms of a set of boson fields. 

In the previous section we showed that the fermion field ip a {x) can be expressed as a 
coherent state of boson operators. The generalization of this boson representation to an 
arbitrary state of the Fock space is surprisingly straight forward. 

As the boson excitations form a complete basis of any N particle subspace we can write 
an arbitrary state as some set of boson excitations above an N particle ground state, 
that is, |*> = /({6tJ)|jV ). 

Using the commutation relations for ip a (x) and , and the relation, 



e 



^ {x) f{{blA)^ {x) = f{{hl a , + <W7^}) (12) 



2vr • , 



Lq 



we can derive the bosonization identity, 



i/) ff (x)\V) = F a \ a {x)e^ {x) e-^ {x) \y) . (13) 

The boson representation of the fermion fields allows us to obtain fermion expectation val- 
ues in terms of the expectation values of the exponentials of the boson fields. In similar 
treatments, (Emery 1979; Mahan 1981; Delft and Scholler 1998), the boson operators b qa 
and b^ qa , are given factors of — % and i, respectively, which changes the exponentials in Eq. 
(|13|) from ± to — i for both terms. 
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B. Bosonization of Hamiltonians 



The definition of the boson operators as bilinear combinations of fermion operators would 
suggest that the two particle interaction terms of a fermionic Hamiltonian, which may be 
expressed as an effective density- density interaction, may be rewritten in terms of a sum of 
quadratic combinations of boson operators. 

If the kinetic term of the fermion Hamiltonian, which in the case of a Fermi gas, only 
contains a single bilinear combination of fermion operators, was also expressible in terms of 
quadratic combinations of boson operators, then the model may easily be diagonalised in 
terms of the boson operators. 

By commuting the boson annihilation operator b qa with the kinetic term of a fermion 
Hamiltonian of the form, H = J2k,<r e (k) ^ka^kai 

[Ho, bj = J2« k ~ ~ < k )) cl- q ,aCk* , (14) 

k 

we find, if the quantity e(k — q) — e{k) is independent of k, the commutation relation is that 
of a lowering operator b qa . The the kinetic term of the Hamiltonian may then be expressed 
as a bilinear combination of the operator b qa and its adjoint. This identification of the boson 
and fermion Hamiltonians is known as Kronig's identity (Kronig 1935). 

By examination we can see, for the term e{k — q) — e(k) to be independent of k, the 
spectrum must be linear in k, that is, e(k) = Ak + B, for some constants A and B. For this 
reason, in the next chapter, we begin our examination of fermions in one dimension with a 
gas of fermions with an infinite linear dispersion. 

III. FERMIONS WITH LINEAR DISPERSION 

These fermions will have an infinite linear dispersion e{k) = vp{k — kp), for vp the Fermi 
velocity, and kp the Fermi momentum. The model has an infinite Fermi sea of positron 
states. However, these states do not contribute to the low energy physics of the model, and 
therefore the physics is believed to be descriptive of the low energy sector of more realistic 
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one dimensional models. Restriction of the particles to movement in one direction also has 
physical realisations. The transport in quantum Hall edge states is an example of one such 
occurrence. 

The construction of a model with infinite linear dispersion provides an exact bosonization 
of the model with a non-interacting boson Hamiltonian. All the properties of the fermionic 
system may then be obtained by direct calculation, using the boson representation. As this 
model is critical at zero temperature, the scale and Lorentz invariance of the system allows 
it to be solved using conformal field theory (Christe and Henkel 1993). This is possible due 
to the Fermi velocity acting as a Lorentz invariant velocity equivalent to the speed of light 
in conformal field theory models. Models which contain excitations with more than a single 
velocity may still be solved using a product of Virasoro algebras and a dressed charge matrix 
(Frahm and Korepin 1990). 



In the non-interacting case, H = vf Sfcl-oo ~~ ^f) cjxfc> the direct solution of the 
fermion model requires less work than the bosonization method. However, the bosonization 
gives some insight into the structure of one dimensional fermion systems. The bosonization 
solution also extends easily to the interacting case, unlike the direct fermionic solution. 
Being on a lattice we also have: 



rik G Z, 5b G [0, 2) a parameter that determines the boundary conditions, and L the length 
of the system. 

Defining the bosons, for q > 0, 



A. The Non-Interacting Spinless Model 




(15) 




(16) 



we then commute the boson operators with the Hamiltonian, to find, 
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[^o, 6,] = Jj% v F ^k> (k' -q- k F )c\, c k , - (k' - k F )c{,_ c k , 



W 

= \ff-VF E fe '(-?)4'- ? % = -v F q b q . (17) 

Using Kronig's identity (Kronig 1935) we can rewrite the Hamiltonian as a sum over bi-linear 
boson operators, 

H = v F J2qbih + C , (18) 
q 

where C is an operator that depends on the energies of the bosonic ground states. To 
determine the operator C we apply the fermion Hamiltonian to an arbitrary ground state 
of N particles, 

H \N ) = b B \N ) = ^N(N + 1 - S B )\N ) . (19) 

To obtain this result we must include in the bosonized Hamiltonian the same combination of 
fermion number operators, which must also be normal ordered, as N does not include states 
below the zero particle ground state energy. With the zero particle ground state denoted 
by |0), and defined by c k \0) = for k > 0, and c k \0) = for k < 0, we define the fermion 
number operators by, 

N = E [c& - <0|4c,|0>] . (20) 

k 

Hence, we may write the full boson Hamiltonian as, 

Ho = v F J2Q KK + ^N(N +1-5 B ). (21) 

q L 

This Hamiltonian describes a system of free bosons with energies given by e q = v F q for 
q > 0, with an additional term that gives the bosonic ground state energies for each fixed 
particle number subspace of the fermion Hilbert space. 



1. Distribution Function for the Bosons 



In the finite temperature calculation of expectation values, we must average over an 
ensemble of states. We must therefore check that the extra term in the Hamiltonian, which 

14 



gives the energies of the "ground states" for the bosons, does not change the distribution 
function of these bosons at finite temperature. 

To begin with we simply calculate the partition function, Z, where, 

Z = En E W (iV, {n q }\e~^\N, {n q }) 
= En e H»=3W+i-'*) £ w (iV, {n q }\e-^^\N, {n q }) 

= En e -/»=?^+i-fe) / n 9 (l - e"^) , (22) 

the third line is obtained from the second, by first expressing the sum in the exponential as 
a product of exponentials, and secondly summing over the occupation numbers. The series 
is summed by using E^Lo x ™ = (1 ~~ a;) -1 . Further, we note that the chemical potential for 
the bosons vanishes. 

Calculating the boson distribution function, we obtain, 

(b\b q ) $ = ±E N E { n q , } (N, {n q/ }\e-^b\b q \N, {n q ,}) 
= \ En e-^^+i-fe) E W} (iV, M\e-^ e < n <b\b q \N, {n q ,}) 
= i En e-^ N ^- s ^ E {v} 11, e-^'V (N, {n q ,}\b\b q \N, {n q ,}) 
= En= ne~^ n (l - e"*«) E W } ^e^W (l - e^v) = l/( e ^ - 1) . (23) 

Above we use the series summation E^Lo nx ™ = ~ x )~ 2 -> an d n °te that the sum over 
states of the product, in the second last line, is equal to one. 

The terms introduced by the fermion number operator in the Hamiltonian cancel, leaving 
a free boson distribution function. The lack of dependence of the distribution function on 
the boson ground states can be understood physically, in that, each particle- hole excitation 
only depends on the difference in energy between the hole and particle, and not on the 
energies of the particle and hole individually. 

The partition function, however, relates to the underlying fermion system, and hence is 
not equal to the partition function of a gas of bosons, and will give fermionic thermodynamic 
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properties. The equality of the fermion and boson partition functions is shown explicitly 
in the article by Haldane (1981), and is essential in proving the completeness of the boson 
states. 



2. Time Dependence of the Boson Fields 

For the non-interacting case, the time dependence of the boson operators is trivially 
obtained due to the non-interacting nature of the boson Hamiltonian. In deriving the time 
dependence of the boson fields we shall generalize the result for the case of a free boson 
Hamiltonian with arbitrary dispersion. This will allow the use of these results in later 
sections. 

The boson operators b q have a time dependence given by, 

b q (t) = e l H q e~ lflt = e- lt " l b q , (24) 

where e q is the dispersion, which in the non-interacting case is simply e q = vpq. 
We define the boson fields, 

<P(?) = ~ E \ T-/ qX \ > <P\*) = - E \ T- n e - iqXb l > (25) 

which when commuted with the free boson Hamiltonian, introduce a factor given by the 
dispersion e q for each b qa operator in the sum, 



— / In - I lir 

[H Mx)] =-E^T-l^ X HW = T.^\ J-/ qX \ ■ (26) 

q,q' V L( l q V L 1 

Thus, the time dependence of the field operators is easily obtained by moving the time 
dependent operators inside the sum, to obtain, 

<p(x, t) = e iAt <p(x)e- iAt = -J2 \h^e iqx e- le « l b q , (27) 

q V Lq 

which in the non-interacting case reduces to (p(x,t) = (p(x — vpt). Similarly, we find, 

lOjr 

V \x,t) = Y.\j Tq e- iqx e^ t b\ (28) 
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and, in the non-interacting case (p^(x,t) = (p^(x — v F t). 

Returning explicitly to the case e q = v F q, the simple time dependence of the <p(x,t) 
boson fields may be shown to extend to the time dependence of the coherent state e~^ x '^ = 
e -ip{x-v F t)_ Thus, the evolution in time of the exponentials of the boson fields results in a 
translation of the fields by v F t. 

Using a different derivation we can show explicitly how the Hamiltonian acts on the 
boson field state as a generator of translations. To begin, the commutators of the fields with 
the Hamiltonian may be rewritten in the form, 

[H , <p{x)\ = iv F ^- , [H , v \x)} = ivF 9 -^- • (29) 

From this we can calculate the commutator of the Hamiltonian with the exponential of the 
(p(x) fields, from which we obtain, 

[H Q ,e-^\ = e-^) (-[H ,tp(x)\ + ±\[Ho,ip{x)],ip(x)\ + ...) 

= - e -**) [H M*)\ = -iv F e-*M2@?- = iv F fe-*M , (30) 

where the second line follows from the first because the boson field commutes with its 
derivatives. Applying the Baker-Hausdorff identity to the time evolved exponential of the 
boson field, 

e iH t e - v ( x ) e -m t = e -„(*) _ Vp vz + / f) 2 e + = e ^ {x - VFt) (31) 

ox 2! ox 

Similarly, we can show, 

[Ho, e* Hx) ] = ivp^-e^W , e iAot e^ {x) e~ iflot = e vHx ' VFt) . (32) 
ox 

The dependence of the fields on the combination x ± v F t, may also be determined from 
a Lorentz invariance of the Hamiltonian. The velocity v F replaces the velocity of light as 
the fundamental velocity in the theory (Senechal 1999). 

The Lorentz invariance remains even in the presence of delta function interactions, due 
to the separation and independence of the charge and spin sectors of the interacting model, 
and linear spectrum of the excitations. In this case, there exist two independent Lorentz 
invariant systems with velocities v e and v x . 
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3. The Klein Factors 



Here we calculate the commutators of the Klein factors F, F\ with the phase factor 
X(x), and also the time dependence of the Klein factors with the free boson Hamiltonian. 

First, we calculate the commutation relations between the Klein factors and the phase 
factor, 

[F, X(x)] = [F, ^fe-^(*-5 fe ) 4 ] = y/^eW^fre- 43 *"*] 
j e i%6 B x e -i%Nx fi^F^ N]-± i^f-x) 2 [[F, N},N} + ... 

= ^eM^>F, (33) 
which gives the relation, 

FX(x) = (l + e-^ x ) X(x)F . (34) 

The time dependence of the Klein factors depends only on the number operator part of 
the free boson Hamiltonian. Thus, 

F(t) = e 'Tf^(W+l-fe)t p e -i^N(N+l-5 B )t 



e 



T~\N+i--)tp = fZ^i^tx^p . (35) 



These phase factors have all previously been ignored in calculating expectation values of 
fermion operators (Mattis and Lieb 1965; Voit 1995; Delft and Scholler 1998), when taking 
the limit L — ► oo. In doing so we are making the presumption that the system has a fixed 
number of particles N, and in the thermodynamic limit N/ L — > p, for some constant p. 

If we take the phase factor for a fixed N, we can see that the phase term 2ttN/ L is just 
the Fermi wave number hp. Extrapolating this to the general contribution of the q = 
phase terms, we can now presume these to generate the chemical potential /i in the boson 
representation. As has been shown in Section [111 A 1|, the explicit chemical potential vanishes 



in the boson representation. Our above reasoning shows that it is retained in the formalism, 
only it is now related to the q = 0, or N operator, dependence. 



To examine the finite size effects on this system, or if the particle number varies over the 
whole fermion Fock space, these terms must explicitly be included. However, in reducing 
a one dimensional system to the linearised model we are only considering excitations close 
to the Fermi surface, which corresponds to obtaining the large distance behaviour of the 
correlation functions and Green's functions. In this case we simply ignore the contribution 
of the Klein operators and phase factors, and add a term e %kpX to account for the chemical 
potential in the correlation functions. 



4- Finite Temperature Green's Functions 

The correlation functions (ip(x, t)^(0, 0))p and (^(0, 0)ip(x, t))p are easily calculated di- 
rectly from the time dependent momentum representation of the fermion operators, followed 
by Fourier transforming over k. From this we obtain, 

The retarded Green's function, Gr(x, t) = —i9(t){[ip(x, i),ip*(0, 0)]+)p, can be calculated 
easily for the free fermion Hamiltonian, and for inverse temperature f3. From this calculation 
we obtain, 

l^ R {X,i)— l^lj L Li [ 1+e -?(t t -(.) l +e /3(e fc -M) 

= -i9(t)^E k e ite e-^ efc -^* = -i9(t)e ikFX 5(x - v F t) . (37) 

The calculation using bosonization is more complicated for this simple case, but the 
generalization to the interacting case is straight forward, whereas the fermionic case is not. 
By expressing the fermion operators ip(x, t) in terms of the boson operators, we note that the 
Klein factors pick up a time dependence due to the number operators in the Hamiltonian, 
and the exponential of the boson operators obtains a time dependence from the bilinear 
boson term. That is, 

ij(x, t) = e ikt i>{x)e- ikt = J^e-'^XivF^FXix^^-^e^-^ . (38) 
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The correlation functions are now given by, 

(i/>{x, t)^{0, 0)){f = (e~ i2 £ VFt ^X{v F t)F\{x)X^{0)F^e {fi ^ x - VFt h-^ x - v ^e-^^e^) i3 

= (l + e~ t2 r x ) e~ i2 r VFt (\{x + c^^^e-^e^'M )) . (39) 

The term involving the fermion number operator will introduce phase factors in the sum over 
the boson ground states, which will not necessarily cancel with the terms in the partition 
function, giving phase factors related to the chemical potential. We ignore these for now, 
and concentrate on the expectation of the boson operators. 

As the time dependence may be included by the exchange of x for x — vpt, from this 
point we only explicitly write the x dependence. 

The expectation of the four exponential boson fields can be simplified by commuting 
the exponentials to move all the boson creation operators to the left of the annihilation 
operators. The commutation introduces a term, which we denote by, 

D(x, 0) = eW*).^)] = exp f ^ ]T —) , (40) 

that diverges in the limit of large system size. Therefore, this term must be included back 
into the expectation at a future stage of the calculation in order to obtain correct results. 

The expectation of the exponentials of the boson fields can be solved simply by using a 
relation for expectation values of exponentials, where, 

^exp (]T A q b^ exp ^ Bjb?J ^ = exp (]T A q B q (e^ - l)^ , (41) 

in the case where the Hamiltonian is identical to that of a simple harmonic oscillator, so we 
can follow it's solution (Mermin 1966). Ignoring the preceding phase factors generated by 
the Klein and A terms, we may express the expectation of the boson fields in the form, 

(i>(xW(P)) P ^ D(x, 0)^ exp (E 9 ^/ff (1 - e-^)b\) 

xexp(E ?v /f(e i9a: -l)6 3 )) . (42) 

which gives the A q and B q terms in Eq. (flU). 
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In the continuum limit the sum in the exponential in Eq. (£|l]) becomes an integral, and 
the correlation functions can be evaluated explicitly by solving the integrals, 



cosqx 



-dq + 



dq 



,iqx 



(43) 



/„ q( e ' 3v Fi - 1) Jo q 

where the first integral gives the contribution of the expectation of the fields, and the second 
integral comes from the D(x, 0) term generated by the commutator. 

By combining the two terms we may be solve the integral explicitly, using the integral 
representation of the logarithm of the gamma function given by Ryzhik and Gradstein 
(1965), and we obtain the expression, 



oo 2 
o q 



cos qx — 1 e 
+ 



tqx 



dq = In 



e Pv F q _l 2 

where the divergent constant Q, is given by, 

Q 



Til 



IX 



j3v F/ 



+ Q 



(44) 



-0v F q 



q 



dq . 



(45) 



The divergent constant Q does not depend on x (or t), unlike the term D(x,t). Thus, it 
may now be treated as a normalization factor that can be essentially ignored, as it plays no 
part in the dynamics or spatial dependence of the correlation functions. 

The exponential of the right hand side of Eq. (fWj) may be rewritten in the form (Spain 
and Smith 1970), 



(46) 



^Fsinh^ 

The above equation again gives us the fermion correlation function, except for the phase 
factor e tkpX contributed by the chemical potential, which we have ignored by omitting the 
leading phase factors in going from Eq. (^) to ([42j). 

The commuting of the two fermion operators within the expectation changes the diver- 
gent factor D(x,t) to D(—x, —t), which merely results in a change in sign of x (and t) in 
the integral of Eq. (f44|). From this the relation, 



(if>\0)if>(x,t))f,= (il>(-x, -t)^(0)) p 



(47) 
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follows immediately. This simplifies the derivation of the retarded Green's functions, which 
contain a sum of these terms. 

The time dependent correlation functions, as well as all the finite temperature Green's 
functions, may now be calculated. To obtain the time dependence of any of these functions 
we simply replace the variable x by x — vpt. 

By introducing a convergence term, e~ aq , the integral may also be solved in closed form 
as a ratio of gamma functions and a term ij (x + ia), which reduce to the form, ir/ sinh7rx, in 
the limit a — > (Emery 1979). This result uses non-normal ordered terms as boson operators 
which require the introduction of a term a^ 1 to retain the fermion field anti-commutation 
relations. The convergence term introduced in the solution of the integral then cancels with 
the a -1 term to produce the desired result. 

Preliminary results in obtaining Eq. (??) for a non-interacting spinless model, via 
bosonization has been reported previously by us (Bowen, Gulacsi and Rosengren 2000). 
However, Eq. ([P]) was known from conformal field theory (Cardy 1984) and it is widely 
recognized (Christe and Henkel 1993) that it is one of the most important results of the 
application of conformal invariance to critical phenomena. 

B. The Interacting Spinless Model 

The interaction term for spinless fermions in one dimensional is of the form, 

Hint = yEE V (l) d k+/k'~g d k' d k , ( 48 ) 

q^O k,k' 

the q = term can be neglected due to the background Coulomb interactions with the 
lattice (Mahan 1981). In terms of boson operators this becomes, 

H INT = WqV(q)(blb q + b q bl) . (49) 

We commute the second term and renormalize to remove the constant interaction energy, 
then we can express the interaction as, 
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H INT = Y,qV(q)b\b q , (50) 

which when added to the non-interacting boson Hamiltonian changes the dispersion of the 
bosons, or if V(q) — V is a constant, the interaction renormalizes the Fermi velocity, 

H = v F J2qbi\ — H Q + H INT = Y^{v F + V{q))qb\b q . (51) 



1. Finite Temperature Correlation Functions 

The same argument used for the non-interacting correlation functions extends to the 
interacting case. The only difference is in the time dependence of the boson operators b q 
and b q , which now have a different dispersion e q in the time dependence. 

The equal time correlation functions are obtained by solving the integral, 

dq . (52) 

where e q is the new dispersion function of the bosons. 

For V(q) = V, the real space interaction is a delta function interaction S(x), meaning 
two particles only interact if they collide with each other. For spinless fermions this is not 
possible, as the Pauli exclusion forbids identical fermions from being in the same state. 
Hence, for a delta function potential in this model the result should be the same as for a 
non-interacting system. 

In the above calculations we see that the delta function interaction changes the finite 
temperature correlation functions by rescaling the Fermi velocity. This is the equivalent of 
rescaling the inverse temperature, with the scaling determined by the coefficient, or strength, 
of the delta function coupling. Apart from the rescaling, however, the distribution function 
n(k,/3), is unchanged from the non-interacting case. 



/ 

Jo 



oo 2 

Q 



cos qx — 1 e iqx 
- 1 + ~2~ 
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C. Interacting Spin-^ Fermions 

Introducing spin to the fermion model fundamentally changes the treatment of the in- 
teraction terms in terms of boson excitations. Maintaining full SU (2) spin invariance in the 
Hamiltonian requires the use of non-Abelian bosonization techniques (Witten 1984). 

In this section, we shall assume coupling is between only the z components of the spin. 
Hence, we add the z spin component to each fermion, which may be in the states denoted 
by | or | . The Hamiltonian is then, 

H = v F J2 J2( k - foOcfeA, + 7EEE v ^'(q)c{ +q ,A'- q yC k , a ,c ka , (53) 

o- k n cr,er' g^O k,k' 

where a, a' G {T, !}■ 

Writing the Hamiltonian in terms of the boson operators b qa , where er now denotes the 
spin index, and V\\(q) and V±(q) denote the parallel spin and anti-parallel spin components 
of the interaction coupling, we find, 

h = E E («*■ + Mi (?)) q hX + E^i (?) {KAi + KM ■ ( 54 ) 

a 1 <J>0 

The parallel and anti-parallel interaction couplings V\\(q) and V±(q) may not be equal, as the 
parallel component must take into account an exchange term due to the indistinguishable 
nature of the scattering fermions. For the anti-parallel case, the scattering fermions are 
now distinguishable and hence no exchange term is present. This may lead to the different 
effective couplings under the bosonization of the Hamiltonian (Mahan 1981). 

1. Spin and Charge Excitations 

Overhauser (1965) has shown that the excitation spectrum of this model is completely 
described by two types of excitations, density oscillations and spin waves. The operators 
that describe these excitations are given by, 

& = 73 (^ + hi) e\ = 73 Kt + Ki) 
X« = 71 {Kl - hi) Xj = 71 {hi ~ hi) > (55) 
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where g q denotes the density excitations, and Xq the spin excitations. Because the particles 
are presumed to be charged, the density of particles also corresponds to the charge density, so 
we may call the density excitations charge excitations. The commutators between different 
operators all vanish, meaning the spin and charge excitations are independent of each other. 
The commutators of the charge and spin operators with the Hamiltonian give the results, 

[H, Q q \ = -(y F + V\\ (q) + Vj_(g)) q g q = -v e (q)q Q q 
[H, Xq] = ~(vf + V\\ (q) - V ± (q)) q % = ~v x {q)q % , (56) 

which again entails we can write the Hamiltonian in the form, 

H = Z v ei ( l) ( l & + Z v x(o)Q x\x q + C , (57) 
q q 

the two types of excitations thus exist independently of each other, and may obey different 
dispersion laws. Any difference in the dispersion means that the two types of excitation may 
separate with time, leading to a phenomena known as spin-charge separation. 

The operator C contains the q = components which account for the fermion number 
operators which may also be split into spin and charge operators. Again, this term relates 
to the ground state energies of the boson states, and represents the chemical potential of 
the fermion system. 

This can be achieved by establishing charge and spin ground states, 

A>^I(A> T + A>J and N x = i(JV T - , (58) 

respectively. These operators may then be separated into the respective spin and charge 
components of the Hamiltonians (Voit 1995). 

2. Bosonization in Terms of Spin and Charge Excitations 

Because the spin and charge operators are constructed from the boson operators b qa , we 
can rewrite the bosonization identity between fermions and bosons in terms of the spin and 
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and hence, 



charge operators. This allows us to easily separate the expectation values into a product of 
spin and charge components. 

We note that, for s = 1, — 1 where a =f, |, respectively, we have, 

\c = -j= (pq + S Xq) , K° = ^ {p\ + S 4) ( 59 ) 

<ft{x) = ~ E g y/^e-** (ej + sx\) = <p),(x) + s^ x (x) . (60) 

As the spin and charge operators commute, we can split any exponentials of the sum of 
these operators into a product of exponentials of each operator. These exponentials will 
then commute past each other, and the single fermion operators may be represented as, 

Mx) = F a \ a (x)e^ (x) e-^ x) e slp * {x) e- slp x {x) . (61) 

The expectation values of the fermion operators then become a product of two boson 
expectation values, with the charge operators traced over a set of charge excitation states, the 
spin operators traced over a set of spin excitation states. The boson excitation states can be 
factored into separate spin and charge number states, as these excitations are independent. 

For single particle correlation functions involving fermions of opposite spin, the Klein 
factors will no longer vanish. The correlations of the most relevance in these models, such 
as singlet and triplet superconductive pairing, and magnetic susceptibility, all involve two 
particle correlation functions in which the Klein factors will again cancel. 

Because the spin-charge boson Hamiltonian is non-interacting even when the underlying 
fermions are interacting, the expectation values obtained are identical to the non-interacting 



fermion case, with only the factor ^J2n/Lq changed to ^Jn/Lq. The dependence on the spin 
state a given by s = ±1 vanishes in the correlation functions, as only factors of s 2 = 1 appear. 
In taking the continuum limit, where the sums in the exponentials become integrals, the 
factor of L/2n does not cancel the ir/L, and each correlation retains an extra power of 1/2. 
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In the case of a real space delta function interaction, we use the Pauli exclusion principle 
to set Vjj = 0, and we have only a delta function interaction V — V±, for fermions of opposite 
spin. The time dependent correlation functions may then be written explicitly as, 

(M*,t)1>i(0,0))e = k/^sinh ^"^ ] 172 [vr/^sinh^p^] 1 ^ , (62 ) 
where v e = vp + V and v x — vp — V. 

IV. THE TOMONAGA-LUTTINGER MODEL 

The Tomonaga-Luttinger model is an exactly soluble model of a one dimensional in- 
teracting electron gas. Tomonaga (1950) developed a model of interacting fermions in one 
dimension in which the density operators had approximately bosonic commutation relations. 
This model contains fermions with linear dispersion and a finite bandwidth. 

In order to make the bosonic commutation relations of the density operators exact, 
Luttinger (1960) extended the model to contain two independent branches of fermions with 
infinite linear dispersion. One branch has a spectrum given by = Uvf and the other by 
efc = —kvF- The two branches are called left movers and right movers, respectively. 

Interactions are restricted to small momentum transfer forward scattering processes in 
order to make an exact solution possible. The effects of adding large momentum transfer 
interactions to the Tomonaga-Luttinger model is the subject of the next chapter. 

The model was correctly solved by Mattis and Lieb (1965) by diagonalisation of the 
Hamiltonian. The zero temperature single particle Green's functions were then calculated 
explicitly by Theumann (1967), and Dover (1968). The representation of the fermion op- 
erators in terms of boson operators, discovered independently by both Luther and Peschel 
(1974), and Mattis (1974), made further calculation of correlation functions straight for- 
ward. The single particle spectral functions were calculated independently by Meden and 
Schonhammer (1992) and Voit (1993). 

The importance of the Tomonaga-Luttinger model lies with the conjecture by Haldane 
(1981) that the model is the asymptotically exact model for gap-less one dimensional inter- 
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acting fermions. The essential physics of each of the gap-less degrees of freedom is then only 
dependent on the values of two parameters, which may be determined by methods such as 
perturbative renormalization or numerical solutions. 

In this chapter we review the bosonization solution for the Tomonaga-Luttinger model. 
We obtain the finite temperature single particle correlation functions. Higher order correla- 
tion functions and Green's functions may be determined by similar methods, and their zero 
temperature form is covered elsewhere (Mattis and Lieb 1965; Voit 1995). The calculation 
of the retarded single particle Green's function at finite temperature for a delta function in- 
teraction is also contained in the article by Emery (1979) as well as some finite temperature 
susceptibilities. 



A. The Luttinger Hamiltonian 

The model contains two species of fermions, left movers and right movers, that are 
represented by the fermion operators c ka , d ka , and their adjoints. The Hamiltonian consists 
of the kinetic term, 

H = v F J2(k - k F ) 4, ff c fcjtr + vJ2(~k - kp) d\ a d k a , (63) 

fc,<7 k,a 

and interaction terms, 

Hint = f J2ay J2 q 92,aa> P c ,AQ)Pd,a'(—Q) 

+ i:J2ay E ? >o#w [pc,*(q)Pc,A-q) + PdA-tiPdAti] ■ ( 64 ) 

The interaction Hamiltonian contains only forward scattering terms, backward and Umk- 
lapp scattering are not included in the Tomonaga-Luttinger model. The coupling coefficients 
are labelled according to the g-ology prescription. The model is a generalization of the 
Tomonaga-Luttinger model, and represents the continuum limit of possible lattice models. 
The kinetic term is equivalent to that of the Luttinger model, with two linear branches of 
fermions. 
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The interactions consist of four terms, with coupling constants labelled by gi, for i = 
1,2,3,4. Each coupling constant may also be spin dependent g iaa i = + gi±, where || and 
_L, refer to parallel and anti-parallel spins, respectively. 

The terms #2 and g^ correspond to forward scattering, or small momentum transfer 
scattering. The g^ term corresponds to two particles on the same branch, and #2 to one 
particle on each branch scattering. 

The backward scattering term g\ relates to large momentum transfer of particles onto 
opposite branches. In models including spin, the parallel component of the backward scatter- 
ing term may be viewed as an effective forward scattering term, as the interaction represents 
the exchange term of the #2 term. In this case the (?2|| term may be rewritten to include 
the contribution. For anti-parallel spins the scattering electrons are distinguishable, and 
hence do not have an exchange term. 

The Umklapp term g% is generally only significant in lattice models when the band filling 
is commensurate with the lattice, that is 4kp is equal to a reciprocal lattice vector. Two 
particles on the same branch may then scatter onto the opposite branch, and this term 
becomes important. 

The density operators in Eq. (Q) have been defined as, 

PcAl) = Z d k+q,aCk,a , Pd,M) = Z ^k+q,Jk,a ■ ( 65 ) 
k k 

The commutators of the density operators are obtained from, 

\PcM)iPc,A-<i]\ = s ™' Z (4+ 9 ,<A +(? v - 4+ a -5',<A<r) • ( 66 ) 

k 

Since the operators must be normal ordered in order to subtract them properly, we find the 
commutator comes out as, 

\PcA-J)iPc,<t(o)] = <W<Wtt ' ( 67 ) 
which gives an approximately bosonic commutation relation. 
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Thus, we can define boson operators to make the commutation relations exact. Even 
though there are now two species of fermions, we can in fact, assign a single species of boson 
operator for all q (Emery 1979). This is done by assigning, for q > 0, 



K = yl^M-s) , b\a = \jj^pc,M > ( 68 ) 



and similarly, 



h - qu = {rq^'^ ' b ~ qa = yl?^ - ^ ■ ^ 

These operators now obey exact boson commutation relations, [b qa , b q/fT ,] = S qq 'S aa '. 

All other boson fields generated as functions of b qa and b qa , in previous chapters may 
now be defined analogously in terms of these operators. It should be noted, however, that 
the boson operators for q > correspond to the right moving fermions, and the q < bosons 
to the left moving fermions, so the boson field operators corresponding to the left and right 
movers must have sums over the appropriate range of momenta q. 

The commutator of the kinetic term of the Hamiltonian with any boson operator gives, 

[ H o, b qa ] = -v F \q\b qa = -e{q)b qa , (70) 

where the dispersion is e(g) is Vpq for q > 0, and — v F q for q < 0. 

From Kronig's relation, the kinetic term can be written in a quadratic form of the density 
operators, 

A) = E E VFQ ib\ a b qa + blj_ q „) ■ (71) 

o- q>0 

Also, because the bosons are defined as excitations above an iV particle ground state, the 
Hamiltonian must include terms that include the energy of the different bosonic ground 
states. These terms are not required for the calculations in this chapter, and are hence 
omitted. 

The interaction term in terms of the boson operators is obtained by direct substitution, 

Hint = T, a ,a' E 9 >o ?2, ro '9 (J>lJ>- qer ' + b qa b_ qa ?) 

(tt qa b qa , + bl qa b_ q(7 >) ■ (72) 
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B. Diagonalising the Hamiltonian 



The eigenvalues and eigenstates of the Luttinger model were first solved by diagonalisa- 
tion of the spinless Hamiltonian by Mattis and Lieb (1965). In this section we diagonalise 
the bosonic Hamiltonian by a Bogolubov transformation (Emery 1979, Mahan 1981, Voit 
1995). 

For each degree of freedom in the model, we obtain a renormalized velocity parameter, 
v, and a renormalized coupling parameter, K, that completely describe the properties of the 
system. 

1. Diagonalisation for Spinless Fermions 

We start by diagonalising the case with spinless fermions, as this generalizes easily to 
the case with spin. The boson Hamiltonian in this case becomes, 

H = J2(v F + g 4 (q))q (b\b q + bl q b_ q ) + ]T g 2 (q)q (b\bl q + b q b_ q ) . (73) 

<3>0 <3>0 

We define the operator, 

S = iT,h {KK + KK) » ( 74 ) 

q>0 

which has the commutation relations [S, b q ] = —i(fi q b_ q and [S, b_ q ] = —i<p q b q , with <p q = <fi_ q . 
The boson operators then transform to give, 

e iSp -is _ cosh <p q _|_ ginh ( f )q 
e i§ b_ q e- i§ = b\ sinh g + b_ q cosh g , (75) 

and the Hamiltonian may be written in diagonal form, 



H D = e iS He~ iS = J2 v f\q\ 
q 

when the condition, 



\ ] 94{q) Y ( 92{qy Z ' 
v F { v F 



blk , (76) 
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is satisfied. This condition may be rewritten in terms of an effective coupling constant K(q), 
which appears in correlation functions, by using the identity, tanh -1 x — \ ln(l +x) — \ ln(l — 
x), to obtain, 

2. The Spin--?; Hamiltonian 

The addition of spin into the Hamiltonian does not change from the spinless case if the 
spins of the bosons are identical, apart from the addition of a sum over spin states. The terms 
that involve boson operators of opposite spins add terms equivalent to those introduced in 
the single branch model considered in Section [III C. 



As done previously, we may obtain the form of the Hamiltonian we need, by transforming 
to the spin and charge variables, 

h = 72 (hi + hi) ' $ = 71 (hi + hi 

% = 72 (hi - hi) > x\ = % ft - hi) ■ ( 79 ) 

Under such a change in variables the interaction couplings transform as, 



92 Q = \ \92\\ + 92±) , 94 e = \ (#411 + 94± / 
92 X = \ (#2|| - 92±) , 9i X = \ (#411 - 94±) ■ (80) 

The transformation then gives a total Hamiltonian consisting of the sum of two Hamiltonians 
of the spinless Luttinger form, as in Eq. ( [T3|) , one for the spin and one for the charge terms. 
These are independent and can be diagonalised by using the sum of the operators, 

S B = 1 Y, (eqQ-q + QqO-q) > 4 = 1 ^x(?) (tfqX-q + XqX-q) • (8 1 ) 

q>0 q>0 
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The diagonal Hamiltonian has spin and charge terms with the renormalized charge and spin 
velocities, 



M = (V + g^(q)) 2 - glM)) 1 ' 2 > Mi) = (V + g^i)? - <4(?)) 1/2 , (82) 



and the renormalized effective coupling constants become, 

^ F + g 4.W- g2g fa) \ 1/2 ( ^ + ^M-^fa) V /2 (83) 
+ 94,(9) + * w ' V«* + »«<(?) + »*(?)/ 



C. Finite Temperature Correlation Functions 

To obtain the single particle correlation functions, we must again express the fermion 
operators in terms of the diagonalised boson operators. The expectation over the fermion 
operators may be rewritten in the form, 

(^ a (x,t)^l(0,0))p = ^Tr e'V^e-^e^fo 0^(0, 0)e~ iS 

= ^Tr e~^ D e i§ ^ a (x, t)^, (0, 0)e~ i§ . (84) 

As the operator e lS is unitary, it can be inserted between operators in the form 1 = e~ lS e lS . 
Also, the trace is invariant under cyclic permutations of the operators. The non-diagonal 
and diagonal Hamiltonians are denoted by H and Hp, respectively. 



1. The Spinless Case 

The correlation functions can be determined by replacing the transformed fermion op- 
erators with the transformed boson field operators. Calculating this for the right moving 
fermion field ipn(x), which is represented by bosons with q > 0, we have, 

x e iS exp (E 9 >o ^¥ q (e iqx ~ 1)6,) e~A , (85) 
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where D(x, 0) is the divergent term generated by ordering the boson fields. The transfor- 
mation leads to terms involving both left and right moving bosons, represented as q > 
and q < terms. These are independent and may be separated, along with the diagonalised 
Hamiltonian, into the product of expectation values for left and right moving fields. The 
right moving fields obtain a cosh0 g dependence, 

( ex P (E - e ~ l9X )H cosh <f> q ) exp (j2 {jTq {etqX ~ 1)h * C ° Sh ^ ) ' (86) 

and the left movers a sinh 0„ term, 



R 



( ex P \E \Z^(! - e~ iqX )K sinh 0, j exp ^ - 1)SL 9 sinh 9 j ^ . (87) 

If the coupling constants are considered independent of the momentum q, which makes 
cf) q = 0, a constant, and the sinh</> and cosh0 terms can be brought outside the summation. 
The form of these correlation functions does not change except for the introduction of an 
extra power term, 1 + 7 = cosh 2 </> or 7 = sinh 2 0, respectively. The term 7 depends on the 
value of the coupling constant K, and is given by, 

7=7^+^-2). (88) 



4 V K 

As can be seen from the definitions, the mixing of fermion species and the anomalous 
power law scalings occur only in the case of scattering between fermions on different branches 
of the dispersion, g 2 7^ 0, which implies if /l. 

With a constant, the divergent term D(x, 0) may be split into two using the identity, 
cosh 2 — sinh 2 = 1, to add a term ±e iqx to the right and left moving components, re- 
spectively. The right moving term is equivalent to the term given in Eq. (0), except for a 
constant 1 + 7, which becomes an exponent. 

The left moving term contains the operators in a reverse order compared to the right 
moving term. This produces an additional term, 2cosgx — 2, which removes the, e tqx , 
component added by splitting the, D(x, 0), term and leaves an x dependent term e~ iqx and 
a constant. Rearranging, we find, 
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oo 2 

o q 



cos qx — 1 e: 



l qx 



dq = In 



e (Sv F q _l 2 

with a new divergent constant Q', given by, 



f3vp J \ {3vf 



+ Q' , (89) 



roc p -Pv F q _ o 

Q' = / dq . (90) 

Ignoring the divergent normalization terms Q and Q', and the phase factor, the equal 
time correlation function may now be written, 

(M*Mo)),* (jjdsTft) ' (91) 

where, a = 1 + 27. This reduces to, a = 1, in the case K = 1. The correlation function for 
left movers t/'l^); is derived analogously, and gives the same result, at least for equal time 
correlations. This result agrees with the conformal field theory results for finite temperature 
(Frahm and Korepin 1990). 

The inclusion of time dependence into the above results requires slightly more care, as 
the terms are obtained from a combination of left and right moving fields. The right moving 
fields maintain the x — Vpt dependence on time, whereas the left moving fields have an 
x + vpt dependence. 

The two properties that determine the non-Fermi liquid behaviour in one dimension 
are spin-charge separation and non-universal power law behaviour (Metzner and Di Castro 
1993; Gulacsi 1997; Voit 1998). We have shown that these two properties are valid also at 
finite temperature. Previously we proved that spin-charge separation is determined only by 
scattering at the same Fermi point. While the above results show that the anomalous power 
law behaviour occurs due to scattering across the two different Fermi points. 

2. The Case With Spin 

To calculate the correlation functions including spin in the model, we must first transform 
into the spin and charge variables, so as to make the Hamiltonian diagonalizable in these 



variables, as shown in Section |I V B 2 . 
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The correlation functions may be expressed as a product of the spin and charge compo- 
nents, where each has the same form as the spinless case presented above. For example, 



'e i§ ^ R (x)^ R (0)e-^) ~ D(x, 0) 

x rw ( e ^ M ex p {^>o v / I( 1 - e ~ iqx )^l) e 

"3>0 A/ Lq 



xe^exp E,> ,& i?1 -1K e-M . (92) 



Under the transformation by the diagonalisation operators, S, the correlation functions 
again split into independent left and right moving components. Each of these components 
is equivalent to a spinless correlation function, with the exception that each correlation 
exponent is now multiplied by a half. The charge and spin components also now have 
separate correlation exponents j e and j x , 

^4(^+7?r 2 ) • ^=H A *+^- 2 ) • < 93) 

and we obtain the result, 



)7m 

for the right movers. 

The time dependence of the correlation functions may be obtained by substituting x+v^t 
and x — Vfjt in each left and right moving term, respectively. Thus the full time dependent 
correlation functions are given by, 

and similarly for the left movers. 



V. THE HUBBARD MODEL 



The Hubbard model was constructed in an attempt to model the correlations of electrons 
occupying the d band of transition metals (Gutzwiller 1963, Hubbard 1963). In contrast to 
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the free electron model of the conduction band electrons, the Hubbard model is a simple 
tight binding model. In a tight binding model each electron is considered as localised on a 
given atomic site, and occasionally the electron will hop to another atomic site, normally a 
nearest neighbour. 

The simplest model assumes each site consists of a single s-band orbital, and may be 
occupied by, at most, one spin up electron and one spin down electron. The interactions 
between the electrons are strongest when they occupy the same atomic site, and all other 
inter-site interactions are ignored. The extended Hubbard model includes a contribution 
from nearest neighbour interactions between electrons. 

As presented in Section |, an exact solution for the Hubbard model been found, where 
the solution depends on an extension of the Bethe ansatz. The ground state energy, wave 
function, and chemical potential for the one dimensional Hubbard model were first found by 
Lieb and Wu (1968). However, the wavefunction is far too complex to obtain expectation 
values and correlation functions by current methods. Takahashi (1972) later formulated the 
thermodynamic Bethe Ansatz equations for finite temperature, from which thermodynamics 
properties of the model may be derived. 

Despite the simple nature of the model the system shows surprising complexity of be- 
haviour. The Hubbard model has also been used extensively in the examination of the 
possible failure of Fermi liquid theory in dimensions higher that one. The higher dimen- 
sional models are not soluble by the same methods as the one dimensional model, and 
techniques such as higher dimensional bosonization are yet to prove useful. 

A. The Hubbard Hamiltonian 

The model has a Hamiltonian of the form, 

H = -J2J2 r Jk (4 A, CT + 4, A*) +^EE + u EE > ( 96 ) 

a j t k cr j a j 

where the operators a^ a and a J cr are the creation and annihilation operators, respectively, 
for an electron in a Wannier state on the atomic lattice site j with spin a = {j, J,}. 
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The first term is the kinetic term, where the numbers Tjk are the transition probabilities 
for an electron to hop from site j to site k. To further simplify the model, the inter-atomic 
hopping can be restricted to nearest neighbour sites, where we write Tjk — r for j — k ± 1, 
and Tjk = otherwise. 

The second term contains the chemical potential //, and the number operator nj >a = 
dj^aj ^, and is used to fix the particle number of the system. 

The final term is the on-site interaction between electrons on the same site. When two 
electrons, of opposite spin, occupy the same site the energy is raised or lowered by an amount 
U, due to forces such as the Coulomb repulsion and exchange forces. The parameter U thus 
gives the interaction strength, and can be repulsive U > 0, or attractive U < 0. 

Transforming the Hubbard Hamiltonian into a momentum representation using plane 
wave modes, 

H + H int = J2(t* ~ It c °s ka) 4,<A,* + u 4+ g , At4- 9 Au > ( 97 ) 

k,a k,k',q 

we obtain the dispersion relation tk = — 2t cos ka, where the hopping matrix element, r, 
determines the bandwidth e max — tmin = 4r. 

B. Bosonization of the Hubbard Model 

Currently, bosonization of the Hubbard model relies on the projection of the model onto 
the Luttinger model. The transformation from a lattice model to the continuum may be 
achieved by a number of methods (Emery 1979, Shankar 1995, Schulz 1993). The Hamilto- 
nian splits into separate bosonic Hamiltonians for the charge and spin sectors. 

Away from half filling the charge sector of the model becomes a Luttinger model. The 
spin sector, on the other hand, contains a Luther-Emery, backward scattering, interaction 
term. For a repulsive on-site interaction, U > 0, this term is irrelevant under renormaliza- 
tion, and does not affect the large distance asymptotic form of the correlations. In the case 
of an attractive on-site interaction, U < 0, the interaction term remains relevant, and a spin 
gap is present. 
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At half filling, the charge sector also gains a Luther-Emery term. In this case the opposite 
applies to the spin sector, a repulsive term, U > 0, generates a charge gap, whereas the 
attractive case, U < 0, is irrelevant under renormalization. The generation of the charge 
gap defines the transition as a Mott type "metal-insulator" transition (Gulacsi and Bedell 
1994a). 



1. Linearisation of the Hamiltonian 

As the methodology of bosonization developed previously does not rely on a continuum 
model, in a first approximation it is possible to directly bosonize the model by linearising 
the dispersion. This gives a Tomonaga model, 

H^v F YX\ k \~ M + U T, E d l+ q ^ k A'- q ,l £ k',i (98) 

k,a q k,k' 

with vp = 2r sin k F a. This may then have positron states appended to give a Luttinger type 
kinetic term. 

The interaction terms away from half-filling are easily obtainable for small positive U, 
where the interactions may be treated as a perturbation. We then have, 

a — U n — U 

92q — 27T > #4(3 — 27T 

92 X = ~£, <?4 X = -g, (99) 

which gives the relations for the exponents, 

/ U \~ 1/2 / U \~ 1/2 , x 

A H 1+ ^) • M 1 -^) • (l00) 

For small values of U/t, these approximations are valid when compared to the exact results 
(Schulz 1993). For strong interactions, or large U/t, the perturbative approximation is no 
longer reliable, and the true values of the correlation exponents may be obtained from the 
Bethe Ansatz solution, or conformal field theory results. 
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2. Physical Fermions 



As we are interested only in the asymptotic low energy properties of the system the 
physical fermions of the model are defined in a way that only includes the components near 
to the Fermi points. We therefore define the physical fermions as, 

4f a (x) ^i>aL(x)+i> aR (x) , (101) 

where the left and right moving fermions are presumed to have a finite momentum cutoff 
around the Fermi points. That is, 

-k F +A 

i> aL (x)~ ]T e ikx d kai (102) 

k=-k F -A 

and similarly for the right movers. 

A representation of the physical fermions in terms of a non-local combination of left and 
right movers is expressed by Voit (1995), where, 

= J + y) + ti R ( x + v)) cot y d v • ( 103 ) 

The solution of the resultant integrals giving the correlation functions obtained from the 
representation in Eq. (|103|) proves to be a difficult task, and the qualitative agreement 
of the simplified approximation with known limits and results does not necessitate their 
calculation, and hence, evaluation is considered a future task. 

C. Correlation and Green's Functions 

In combining operators for the physical fermions we obtain four combinations of left and 
right movers, 

(0) * Vv L (x)^ L (0) + $„r{x)$I l (0) + $ aL (x)$, R {0) + ^ R (x)ti R (0) . (104) 

The terms containing mixtures of left and right moving fermions correspond to excitations 
with q ~ 2fcp, whereas the terms with two of the same fermion type relate to the q ~ 
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components. We can therefore neglect the 2/c F terms as not relevant to the low energy 
physics. 

The asymptotic form of the finite temperature correlation functions for the physical 
fermions as a sum of left and right moving terms, becomes, 



= e 



ikpx ' 



7m 



n ™^s 7 Lh^ W^w-j ■ (105) 

The correlation functions allow us to calculate the single particle Green's function for 
the Hubbard model, from, 



Gl(x,t)^-i9(t)({^ a (x,t),¥ a (0)Y /fi 
= -i6(t) ((* a (x, tM(0))p + -tMm?) , (106) 



we obtain, 



G^ a (x,t) ~ —i9(t) coskFX 



7M 1 + 7M 

2 / \ 2 



n 



7m 1 + 7m 

2 / \ 2 



(107) 



where each term in the sum is in qualitative agreement with the results from conformal field 
theory (Frahm and Korepin 1990). 

In taking the zero temperature limit for the equal time Green's function we recover the 
form, 

G aa (x) ~ x~ l ~ a , (108) 

in this case for a = 1 + 7 e + j x , which is in agreement with the known results (Kawakami 
and Yang 1991, Christe and Henkel 1993, Gulacsi 1997). 
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D. Bosonization and Conformal Field Theory 



The results obtained from conformal field theory are expressed as the correlations of 
the conformal fields themselves, and the physical fields cannot easily be related to these. 
However, the exponents of the leading terms of the asymptotic expansion of the physical 
fields may be approximated by the conformal method. 

The almost identical relation between the correlations of the fermion operators in our 
bosonization approach, and the correlations of the conformal fields may appear to be due 
to the simple relation we have used between the physical fermions and the constructed 
left and right moving fermion fields. The left and right moving fields are independent and 
conformally invariant, each with a separate fundamental velocity. The relationship between 
the physical fermions and left and right moving terms may therefore give insight into the 
relationship of conformal and physical fields in the conformal approach. 

The agreement between the long distance form in the bosonization and conformal ap- 
proaches shows that the bosonization approach is appropriate for obtaining the asymptotic 
correlations. Models that contain gaps in their excitation spectrum are no longer confor- 
mally invariant, and therefore cannot be solved using the methods of conformal field theory. 
The extension of the bosonization method to systems with gapped excitations, such as the 
Luther-Emery model, may therefore give the asymptotic form of the correlation functions to 
a good approximation, where they cannot be obtained from the conformal approach. This 
avenue is currently being explored. 

E. The Metal-Insulator Transition at Finite Temperature 

The metal-insulator transition has been examined at zero temperature (Gulacsi and 
Bedell 1994a), using bosonization of a two band model. The lower band is the Luttinger 
band, whilst the upper band contains the Luther-Emery term. A zero temperature single 
particle Green's function has been proposed for the Hubbard model, by the same authors, 
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again using the two band bosonization method (Gulacsi and Bedell 1994b). 

Given a result for the finite temperature Green's function, it may be compared to the zero 
temperature case proposed by Gulacsi and Bedell (1994a). The temperature dependence of 
the Mott transition in the Hubbard model may then be examined. From the exact Luther- 
Emery model, and solutions in various limits (Voit 1998) has proposed an ansatz form for the 
single particle spectral density function of models in this class. From this proposed function 
it may be possible to calculate the corresponding finite temperature Green's functions to be 
used as an guide. 

Using our bosonization method we have so far been unable to calculate the Green's 
functions of the Luther-Emery model, although the work is still in progress. This is due 
to the complexity of the representation of the spinless fermion operators in terms of the 
physical fermions or their boson representation. 

VI. CONCLUSIONS 

The Luttinger liquid characteristics describe the ground state properties of interacting 
one dimensional systems. All the experimental data, however, are finite temperature mea- 
surements. Because of this it has become timely to comprehensively describe the finite 
temperature characteristics of Luttinger liquids. As a first step in achieving this goal, we 
extended the well-known bosonization method of one dimensional electron systems to finite 
temperatures. 

We have reviewed the derivation of the bosonization operator identity in one dimension. 
The formalism was then used to examine the properties of a simple linear-dispersion fermion 
model, both with and without spin. The existence of spin-charge separation was shown for 
the interacting case, showing that the interactions between fermions at the same point of 
the Fermi surface is the mechanism behind this phenomenon. 

The single particle correlation functions, from which Green's functions may be derived, 
were calculated explicitly to avoid problems of normal ordering of the operators in exponen- 
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tials. Previous calculations have relied on the definition of non-normal ordered boson fields, 
which required a convergence factor to allow solution of the resultant integrals. In using the 
normal ordered representation of von Delft and Scholler (1998), and the iterative solution of 
expectations given by Mermin (1966), we were able to explicitly retain apparent divergences 
and recombine these terms to obtain our results. 

The solution to the Tomonaga-Luttinger model was examined, to show the generation of 
anomalous dimensions of operators which signal the non-universal power law behaviour in 
correlation functions. The dependence of this phenomenon was shown to rely on the scatter- 
ing of fermions of opposite spin on separate branches on the Luttinger model. Correlation 
functions were again explicitly derived, where a mixing of the correlations between separate 
fermion species is observed. 

The asymptotic form of the finite temperature correlation and Green's functions for the 
Hubbard model were calculated. As far as we are aware, these have not previously been 
derived using the bosonization method, and our results agree with the form obtained from 
conformal field theory. It appears, however, that the nature of the relationship between the 
physical fermions of the theory and the constructed boson and fermion fields requires further 
examination. 

In determining the finite temperature behaviour of a two band Hubbard model we hoped 
to model the transition between a Luttinger liquid behaviour at high temperature and a 
Luther-Emery liquid at low temperature. At high temperatures it is assumed that the 
thermal excitations will make the gap in the spectrum of the excitations irrelevant, and the 
model will exhibit Luttinger liquid behaviour, but at lower temperature the gap becomes 
relevant and the behaviour undergoes a phase change to a Luther-Emery liquid. 

Whether the high temperature behaviour of a gapped degree of freedom indeed behaves 
as a Luttinger liquid with perhaps renormalized parameters is the desired goal of this line 
of enquiry. The generalization of finite temperature behaviour of the Luther-Emery model 
to other models within the class is presumed from the renormalization calculations showing 
agreement for this class of models at zero temperature. 
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Difficulties appear with the lack of a simple representation of the original fermion oper- 
ators, in terms of the new spinless fermions, as well as the restriction of refermionization to 
bosonized Hamiltonians with a linear spectrum. 

Further examination of the correlation functions for the Luther-Emery model may show 
whether the iterative solution is tractable in solving the single particle correlator. This 
would allow a simple substitution into the two band model for the Mott metal-insulator 
transition. Due to the equivalence of models, this would also give the single particle finite 
temperature correlation function for the sine-Gordon and massive Thirring models, and 
hence the importance of such solution may be seen. 

The Hubbard model is soluble at finite temperature by use of the Thermodynamic Bethe 
Ansatz. A comparison of the finite temperature properties of the Hubbard model, and those 
of the Tomonaga-Luttinger model, and possibly the Luther-Emery model, may be used to 
clarify the Luttinger liquid conjecture for systems at finite temperature. 
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